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HIGHLIGHTS 


•  Detailed  efficiency  analysis  of  a  shipboard  dc  hybrid  power  system  is  carried  out. 

•  An  optimization  algorithm  is  proposed  to  minimize  fuel  consumption  under  various  loading  conditions. 

•  An  online  optimization  strategy  for  fuel  saving  is  implemented  in  MATLAB/Simulink. 

•  An  offshore  support  vessel  (OSV)  is  simulated  over  different  modes  under  online  optimization  control. 

•  The  method  achieves  around  15%  fuel  saving  with  optimally  utilized  energy  storage  in  the  studied  dc  system. 
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Hybrid  power  systems  with  dc  distribution  are  being  considered  for  commercial  marine  vessels  to 
comply  with  new  stringent  environmental  regulations,  and  to  achieve  higher  fuel  economy.  In  this  paper, 
detailed  efficiency  analysis  of  a  shipboard  dc  hybrid  power  system  is  carried  out.  An  optimization  al¬ 
gorithm  is  proposed  to  minimize  fuel  consumption  under  various  loading  conditions.  The  studied  system 
includes  diesel  engines,  synchronous  generator— rectifier  units,  a  full-bridge  bidirectional  converter,  and 
a  Li-Ion  battery  bank  as  energy  storage.  In  order  to  evaluate  potential  fuel  saving  provided  by  such  a 
system,  an  online  optimization  strategy  for  fuel  consumption  is  implemented.  An  Offshore  Support 
Vessel  (OSV)  is  simulated  over  different  operating  modes  using  the  online  control  strategy.  The  resulted 
consumed  fuel  in  the  simulation  is  compared  to  that  of  a  conventional  ac  power  system,  and  also  a  dc 
power  system  without  energy  storage.  The  results  show  that  while  the  dc  system  without  energy  storage 
provides  noticeable  fuel  saving  compared  to  the  conventional  ac  system,  optimal  utilization  of  the  energy 
storage  in  the  dc  system  results  in  twice  as  much  fuel  saving. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Increasing  concerns  on  environmental  issues  and  fuel  economy 
have  forced  the  maritime  transport  sector  to  pursue  low  emission 
and  fuel  efficient  solutions.  Power  electronics  plays  a  major  role  in 
optimization  of  the  conventional  drivetrains.  In  marine  vessels, 
power  electronic  converters  made  a  breakthrough  with  regard  to 
fuel  saving,  when  they  enabled  electrification  of  propulsion  sys¬ 
tems  through  Variable  Speed  Drives  (VSDs)  in  1990’s  [1].  The  ad¬ 
vantages  gained  from  electric  propulsion  systems  have  ever  since 
promoted  the  concept  of  all-electric  ship,  which  provides  a  com¬ 
mon  electrical  platform  to  supply  the  propulsion  power  and  ship- 
service  loads  [1—3].  Although  the  overwhelming  majority  of  pre¬ 
sent  all-electric  ships  use  ac  distribution  systems,  the  penetration 
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of  power  electronic  converters  into  power  systems  is  giving  rise  to  a 
tendency  towards  dc  distribution  systems  [4-6].  The  incentive  for 
this  transition  lies  on  the  one  hand  in  the  challenges  associated 
with  the  conventional  ship  ac  power  systems,  and  on  the  other 
hand,  in  an  increasing  interest  in  integration  of  energy  sources  and 
storage  devices  with  dc  outputs  [5,7].  The  challenges  of  the  con¬ 
ventional  shipboard  ac  systems  include  the  need  for  synchroniza¬ 
tion  of  the  generation  units,  reactive  power  flow,  inrush  currents  of 
transformers,  harmonic  currents,  and  three-phase  imbalances. 
Fig.  1  shows  single  line  diagram  overview  of  a  shipboard  ac  power 
system  in  comparison  with  that  of  a  shipboard  dc  network.  As  the 
bulky  ac  switchgears  and  transformers  are  removed  in  the  ship¬ 
board  dc  grids,  the  dc  systems  provide  advantages  such  as  space 
and  weight  savings,  and  flexible  arrangement  of  equipment  [8].  Not 
requiring  synchronization  of  generation  units,  a  dc  power  system 
enables  the  prime  movers  to  operate  at  their  optimal  speeds, 
leading  to  significant  fuel  saving.  Fig.  1  also  shows  that  the  dc 
system  eliminates  the  need  for  multiple  stages  of  conversion  and 
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Fig.  1.  Single-line  diagram  overview  of  shipboard  power  systems,  a)  ac  network,  b)  dc  network. 


transformation  that  the  ac  system  usually  requires  when  inte¬ 
grating  a  device  with  dc  output,  such  as  an  energy  storage  device. 

Despite  the  long  history  of  hybridization  and  fuel  saving  prac¬ 
tices  in  land  vehicles  [9-19],  marine  hybrid-electric  vessels  are 
almost  at  the  beginning  of  their  journey,  particularly  with  regard  to 
integration  of  energy  storage  [6,20-28].  Presently,  dc  hybrid  power 
systems  are  being  considered  for  certain  categories  of  ships  that 
experience  frequent  load  variations  and/or  long  time  of  operation 
under  non-rated  conditions.  Offshore  Support  Vessels  (OSVs)  form 
a  major  category  of  interest  [8,24].  The  first  OSVs  with  dc  power 
systems  were  delivered  in  2012,  including  “Jaguar”  with  a  power 
system  from  Alewijnse  [29].  ABB  and  Siemens  will  launch  OSVs 
with  dc  power  systems  in  the  current  year.  “Edda  Ferd”  introduced 
by  Siemens  will  be  the  first  OSV  with  a  dc  hybrid  power  system 
which  benefits  from  variable-speed  diesel  engines  and  batteries, 
and  assures  a  significant  fuel  saving  compared  to  the  conventional 
diesel-electric  systems. 

Although  extensive  research  has  been  reported  on  efficiency 
optimization  of  hybrid  electric  vehicles  [30-33],  this  issue  for  all¬ 
electric  ships  has  largely  remained  unaddressed,  leading  to  the 
research  reported  in  this  article.  This  paper  is  aimed  at  efficiency 
analysis  and  optimization  of  shipboard  power  systems  and  evalu¬ 
ation  of  their  fuel  saving  when  utilizing  the  dc  distribution  and 
incorporating  an  energy  storage  device.  In  Section  2,  fuel  con¬ 
sumption  of  the  prime  mover  is  discussed  and  modeled.  The  effi¬ 
ciency  of  the  electrical  sources,  converters,  and  the  energy  storage 
are  analyzed  and  modeled.  In  order  to  evaluate  potential  fuel  saving 
by  this  system,  in  Section  3,  an  optimization  algorithm  is  proposed 
for  ship  applications  where  multiple  generator  sets  should  be 
optimized  in  terms  of  their  average  powers  and  power  ripples.  The 
optimization  algorithm  extends  the  work  of  [34]  that  suits  HEVs 
with  a  single  generator,  to  all-electric  ships  running  on  multiple 
generator  sets.  The  optimization  algorithm  minimizes  the  fuel 
consumption  by  optimal  power  sharing  among  the  energy  storage 
and  prime  movers.  In  Section  4,  an  online  optimization  strategy 
based  on  the  proposed  algorithm  is  applied  to  an  offshore  support 
vessel  with  seven  different  operational  modes.  MATLAB/Simulink  is 
used  to  simulate  the  system.  The  control  performance  deficiencies 
resulted  from  the  potential  approximations  of  the  online  control  is 
taken  into  account  to  more  realistically  estimate  the  optimized  fuel 
saving  onboard  an  actual  vessel.  The  simulation  results  of  the  on¬ 
line  optimization  control  are  presented  and  compared  with  those  of 
theoretical  optimization  algorithm.  In  Section  5,  the  conclusions  of 


the  study  are  drawn.  The  study  results  show  substantial  fuel  saving 
for  the  optimally  controlled  dc  hybrid  power  system  in  comparison 
with  the  conventional  ac  system  and  also  with  the  dc  system 
without  energy  storage. 

2.  Efficiency  analysis  and  modeling 

Fig.  2  depicts  the  architecture  of  the  studied  dc  power  system  in 
more  detail.  The  power  system  consists  of  four  generation  units,  an 
energy  storage  system  based  on  Li-Ion  batteries,  propulsion  loads, 
ship  service  loads,  and  auxiliary  loads.  In  each  generation  unit,  a 
prime  mover  is  coupled  to  a  wound-rotor  synchronous  generator. 
The  ac  power  of  the  generator  is  converted  to  dc  power  through  a 
three-phase  diode  rectifier.  The  rotational  speed  of  the  prime 
mover  is  adjusted  by  a  local  controller  to  minimize  fuel  consump¬ 
tion.  This  is  not  possible  in  conventional  ship  ac  power  systems 
wherein  the  speed  reference  is  constant  to  maintain  synchroniza¬ 
tion  among  generators.  In  order  to  manage  power  sharing  among 
generators  while  regulating  the  dc  voltage  level,  the  exciters  of 
different  generation  units  can  be  coordinated  by  voltage  droop  with 
secondary  level  control.  The  energy  storage  system  consists  of  a  Li- 
Ion  battery  bank  and  a  bidirectional  dc-dc  converter  that  enables 
to  charge  and  discharge  the  battery  in  a  controlled  manner. 

In  this  section,  the  generation  units  are  analyzed  and  modeled  in 
terms  of  fuel  consumption  of  the  prime  movers,  and  efficiency  of 
the  generator-rectifier  systems.  The  efficiencies  of  the  battery 
bank  and  the  dc-dc  converter  are  also  analyzed  and  modeled.  As 
the  electrical  load  profile  of  the  studied  system  is  available,  the 
efficiency  models  of  the  system  loads  are  not  required. 

2.1.  Prime  mover 

The  main  source  of  power  onboard  an  all-electric  ship  typically 
comprises  a  synchronous  generator,  driven  by  an  internal  com¬ 
bustion  engine  (ICE)  that  is  often  fueled  with  diesel  or  heavy  fuel  oil 
[1  ].  The  specific  fuel  consumption  curves  for  an  ICE,  obtained  from 
engine  efficiency  maps,  have  the  qualitative  behavior  of  Fig.  3  [7]. 
The  combustion  engines  are  continuously  being  improved 
regarding  efficiency  and  emissions;  however,  their  optimum  utili¬ 
zation  factor  of  fuel  is  still  around  40%,  and  the  rest  of  the  stored 
energy  in  the  fuel  is  dissipated  through  heat  or  the  exhaust. 

As  seen  in  Fig.  3,  the  fuel  utilization  of  an  ICE  can  be  improved  by 
adjusting  the  engine  speed  based  on  the  engine  load.  This  implies  a 
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Fig.  2.  Detailed  architecture  of  the  shipboard  dc  power  system  under  study. 


disadvantage  of  the  conventional  shipboard  ac  power  systems, 
wherein  all  generator  sets  must  be  synchronized  at  a  fixed  fre¬ 
quency,  preventing  from  the  optimized  speed  control.  A  dc  power 
system,  however,  enables  the  engines  to  rotate  independently, 
bringing  about  the  possibility  of  optimal  speed  adjustment.  In  ac 
power  systems,  independent  rotation  of  the  engines  could  be  also 
provided  by  use  of  back-to-back  converters  for  each  generator. 
Nevertheless,  this  imposes  more  cost  and  complexity  in  compari¬ 
son  with  the  transition  to  the  dc  system.  In  addition,  considerable 
ship  electrical  power  is  converted  to  dc  to  be  used  by  the  ship 
propulsion  drives,  which  are  the  major  loads  in  many  applications. 
The  use  of  back-to-back  converters  adds  to  the  conversion  stages 
between  ac  and  dc,  which  reduces  the  drivetrain  efficiency. 
Therefore,  the  dc  system  has  become  more  appealing  to  the  ship 
industry  with  regard  to  fuel  saving. 

The  hourly  fuel  consumption  curve  for  an  ICE  can  be  approxi¬ 
mated  by  a  quadratic  function  of  the  load  as  expressed  in  (1).  This 
equation  is  valid  for  either  variable-speed,  or  fixed-speed  opera¬ 
tion,  using  different  sets  of  polynomial  coefficients  [34]. 

Cl  =  C0  +  aPm  +  bPl.  (1) 

where  Co,  a,  and  b  are  the  coefficients  of  the  second-order  poly¬ 
nomial  function,  and  Pm  is  the  engine  mechanical  power.  The 
support  points  for  the  specific  fuel  consumption  of  a  300  kW  diesel 


engine  are  given  in  the  Appendix  [see  l.a)]  under  both  variable  and 
fixed-speed  operations.  The  specific  fuel  consumption  is  the  hourly 
consumption  per  unit  of  power,  which  is  obtained  by 

SFC  =  C0/Pm  +  a  +  bPm.  (2) 

The  above  relations  of  fuel  consumption  also  hold  in  dynamic 
conditions  [34]. 

2.2.  The  synchronous  generator— rectifier  system 

In  a  shipboard  power  generation  system,  the  prime  mover  is 
coupled  to  an  electric  machine,  typically  a  synchronous  generator, 
to  convert  the  mechanical  energy  into  the  electrical  form.  When  dc 
distribution  is  considered,  the  three-phase  ac  power  is  to  be  con¬ 
verted  to  dc  power  through  three-phase  rectifiers.  As  bidirectional 
power  flow  is  not  required  in  such  an  application,  three-phase 
diode  rectifiers  are  appealing  options  due  to  their  structural 
simplicity,  high  energy  efficiency,  and  low  cost.  They  can  readily  be 
used  with  wound-rotor  synchronous  generators  [4,6]. 

Three-phase  rectifiers  are  already  used  in  the  conventional 
electric  propulsion  systems  to  supply  the  dc  links  of  propulsion 
electric  drives.  The  transition  from  ac  to  dc  can  be  viewed  as  the 
enhancement  of  these  dc  links  to  form  a  distributed  dc  bus,  which 
leads  to  elimination  of  bulky  ac  switchboards,  ac  transformers,  and 
their  associated  power  losses  [7].  In  addition,  it  will  be  shown  in  the 
following  that  the  efficiency  of  the  generator  is  improved,  when  the 
prime  mover  is  controlled  at  its  optimal  speed. 

The  main  losses  of  the  synchronous  generator  consist  of  copper 
losses,  iron  losses,  and  mechanical  losses.  The  copper  losses  are 
resulted  from  the  flow  of  current  through  the  stator  and  rotor 
windings,  which  are  calculated  by 

pcu  =  rsI2s  +  ffdJfcj  =  rs(i*d  +  z'sq)  +  rfdl’fd-  (3) 

where  conventional  notations  are  used.  The  iron  losses  are 
composed  of  eddy  and  hysteresis  losses,  which  can  be  obtained  by 
the  following  formula  [35] 

PFe  =  CFe<44  =  CFe^  (Xl dftd  +  *mq'mq)  •  (4) 

where  /3  =  1.5 -1.6.  This  formula  is  represented  in  the  machine 
equivalent  circuit  by  the  core  resistance  Rc  in  parallel  to  the 
magnetizing  reactors  Xmd  and  Xmq.  The  mechanical  losses  are  due 
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Fig.  4.  Overview  of  the  variable-speed  control  of  a  diesel-generator  unit. 


to  friction  and  air  drag  losses,  which  are  proportional  to  the  square 
of  rotational  speed. 

^l.m  Qn^nri-  (5) 

Power  losses  in  diode  rectifiers  are  mainly  resulted  from  the 
conduction  losses,  and  their  switching  losses  are  negligible.  The 
conduction  losses  are  given  by 

Prec  =  2Vfw/dc  + -Pd  onFj2/^..  (6) 

where  VfW  is  the  diode  threshold  voltage,  and  F\  is  the  current 
waveform  factor,  which  can  readily  be  calculated  for  the  typical 
diode  current  waveforms  by  F\  =  /D.rms/fD.avg-  Assuming  that  each 
diode  has  a  half-sine  current  waveform,  we  calculate  7r/2  for  the 
current  waveform  factor. 

In  order  to  evaluate  the  power  losses  in  the  generator-rectifier 
system  and  to  obtain  its  efficiency  curve,  the  simulation  platform 
introduced  in  Ref.  [6]  is  employed.  The  losses  are  investigated  over 
a  wide  range  of  load  for  two  conditions:  1)  when  engine  speed  is 
fixed  at  its  rated  value,  and  2)  when  engine  speed  is  optimally 
adjusted  to  consume  minimum  fuel.  Fig.  4  shows  the  overview  of 
the  variable-speed  control  for  one  diesel  engine.  The  optimal  speed 
against  the  engine  load  can  be  deduced  from  the  engine  manu¬ 
facturer’s  datasheet.  For  constant-speed  operation,  the  reference 
speed  is  fixed  at  its  rated  value.  The  parameters  of  the  diesel- 
generator  sets  are  given  in  the  Appendix  [see  2)]. 

Fig.  5  depicts  power  loss  components  of  the  generator— rectifier 
system  under  fixed  speed  operation,  while  Fig.  6  shows  the  same 
components  when  the  speed  is  varied  to  attain  minimum  SFC.  In 
comparison  of  the  two  figures,  it  is  observed  that  the  generator  core 
and  mechanical  losses  reduce  in  variable-speed  operation,  while 
the  copper  losses  slightly  increase.  Nevertheless,  the  total  loss 


reduction  outweighs  the  increase  in  copper  losses,  which  results  in 
higher  efficiencies  in  variable-speed  operation.  Fig.  7  depicts  the 
efficiency  in  the  two  operation  modes,  wherein  the  difference  of 
the  two  efficiency  curves  increases  as  the  load  decreases. 

The  specific  fuel  consumption  is  calculated  in  terms  of  electrical 
load  by  using  the  efficiency  curves  of  Fig.  7  and  the  SFC  relation  in 
(2).  Fig.  8  depicts  the  curves  of  these  new  SFC  relations  for  fixed  and 
variable-speed  operation,  which  are  fitted  to  quadratic  functions 
with  coefficients  given  in  the  Appendix  [see  l.b)].  These  quadratic 
functions  will  be  used  in  the  development  of  the  optimization 
algorithm. 

2.3.  The  bidirectional  dc—dc  converter 

In  a  dc  hybrid  power  system,  an  energy  storage  device  is 
prominent  to  improve  the  transient  performance  of  the  system,  and 
to  enable  optimal  load  sharing  among  the  sources.  A  bidirectional 
converter  is  normally  used  to  incorporate  the  energy  storage  in  the 
power  system.  Various  topologies  of  bidirectional  dc-dc  power 
converters  are  discussed  and  compared  in  Refs.  [36,37].  Topological 
circuit  of  a  full-bridge  bidirectional  converter  is  depicted  in  Fig.  9. 
The  high  frequency  transformer  is  used  to  provide  high  voltage 
transfer  ratio  to  link  a  low  voltage  battery  to  the  dc  bus.  The 
additional  cost  and  losses  introduced  by  the  high  frequency 
transformer  can  be  justified  when  the  voltage  ratio  between  LV  and 
FIV  is  so  high  that  devices  are  not  economical  to  tolerate  both  high 
voltage  and  high  current  simultaneously  [36].  The  advantage  of  the 
full-bridge  topology  is  to  provide  lower  voltage  and  current 
stresses,  resulting  in  significant  decrease  in  losses  compared  to 
half-bridge  and  push-pull  topologies  [36].  Nevertheless,  it  needs 
high  number  of  switches.  The  power  losses  of  the  converter  of  Fig.  9 
for  its  two  modes  of  operation  are  discussed  in  the  following. 
The  parameters  of  the  studied  converter  are  given  in  the  Appendix 
[see  3)]. 
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Fig.  5.  Power  losses  of  the  generator-rectifier  under  fixed-speed  operation. 


Fig.  6.  Power  losses  of  the  generator- rectifier  under  adjusted-speed  operation. 
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Fig.  7.  Efficiency  curves  of  the  generator-rectifier  under  fixed  and  variable-speed 
operation. 


Fig.  8.  SFC  curves  of  the  diesel-generator  set  in  terms  of  electrical  power  under  fixed 
and  variable-speed  operation. 


The  converter  losses  are  comprised  of  the  conduction  and  the 
switching  losses  of  the  transistors,  the  conduction  losses  of  the 
diodes,  and  the  resistive  losses  of  the  inductor,  the  transformer,  and 
the  capacitors.  The  transformer  core  losses  are  neglected,  as  the 
switching  frequency  is  relatively  low.  In  the  following,  the  con¬ 
verter  losses  in  buck  mode  are  formulated.  The  losses  are  preferably 
derived  in  terms  of  HV-side  power  so  as  to  provide  an  efficiency 
curve  against  exchanged  power  with  the  dc  bus.  Interested  readers 
are  referred  to  Ref.  [38]-Ch.  3  for  details  on  the  derivation. 

The  transistors  conduction  losses  are  obtained  by  (7),  assuming 
that  the  inductor  current  is  ripple-free. 


^Sw,on  —  ^S,orJs,rms  —  ^ 


S,on 


PfwVD 

V, 


hv 


(7) 


Assuming  that  the  transistor  output  capacitance  is  linear,  we  can 
calculate  the  switching  losses  by 


1  Sw.sw 


/sv 


1 


m  4  ^  C0V£m  =  ^/swQjlC 


(8) 


The  following  equations  hold  for  the  diode  conduction  losses 
over  a  switching  period. 


jj  _  (  ^7fwh  +  ,on^L’  (0  <  t  <  DT ) 

PD  ~  \  VfwlL  +  (2i?D,on  ||2RD,on)'L-  iPT<t<  T) 

Taking  an  average  over  the  switching  period,  it  yields 
Pd  =  vfwlL  +  ^d.oi^lD  +  2J?D,on/2(l  -  D) 


(9) 


= 


fw 


f  ft^hv 
\DVhv 


(2  —  D)RDc 


nP i 


hv 


DV} 


hv 


(10) 


The  resistive  losses  of  the  passive  elements  are  obtained  by 


Pi  =  n^L.rms  =  rL 


nP, 


hv 


DV\ 


hv 


Pci  =  rc4ci,rms  =  rCl 


(1  ~D)Phv 


DV, 


hv 


(11) 

(12) 


P  r  1 2  r  (O-Omf2 

PC2  =  rC2/C2,rms  = 


(13) 


Pt  =  rTlhms  =  rT^0=j  .  (14) 

The  power  loss  components  in  buck  mode  are  depicted  in  Fig.  10, 
using  the  above  formulas.  Analysis  of  the  boost  mode  is  carried  out 
in  a  similar  way  to  the  above  formulation,  and  results  in  the  loss 
components  depicted  in  Fig.  11.  For  consistency,  the  formulas  for 
boost  mode  are  summarized  in  the  Appendix  [see  4)]. 

Fig.  12  shows  the  converter  efficiency  curves  for  the  two  modes, 
which  are  slightly  different.  A  simplified  model  of  the  converter 
losses  is  desired  for  the  analysis  of  optimal  conditions  discussed 
later  in  the  paper.  Therefore,  an  averaged  efficiency  curve  between 
the  two  modes  is  considered,  as  the  converter  is  equally  operated  in 
both  modes  to  sustain  the  battery  charge  level.  The  converter  los¬ 
ses,  thereafter,  can  be  modeled  by  an  equivalent  resistance,  as 


Fig.  9.  The  topological  circuit  of  the  full-bridge  bidirectional  dc-dc  converter  under  study. 
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Fig.  10.  The  power  loss  components  of  the  bidirectional  converter  in  buck  mode.  12.  Efficiency  curves  of  the  bidirectional  converter  in  buck  and  boost  modes. 


expressed  by  (15).  This  equation  also  includes  a  small  constant 
term,  resulted  from  constant  switching  losses,  and  capacitor  losses. 
As  seen  in  Fig.  13,  the  loss  model  of  (15)  complies  with  the  detailed 
converter  losses. 

h.bid  =  S^hv  +  ho-  (15) 

Ubus 

ftbid  is  the  equivalent  resistance  of  the  bidirectional  converter.  In 
order  to  provide  a  total  loss  model  of  the  energy  storage  system, 
battery  losses  can  be  integrated  in  the  converter  loss  model  as 
discussed  in  the  following. 

2.4.  The  battery  model 

To  develop  an  optimization  control  strategy,  using  a  reference 
model  of  the  battery  is  essential.  Electrochemical  battery  models 
[39-41]  can  often  provide  high  accuracy  by  capturing  all  key  be¬ 
haviors  of  the  electrochemical  cells.  However,  being  based  on  par¬ 
tial  differential  equations  with  large  number  of  unknown 
parameters,  they  could  lead  to  unnecessary  complex  models  for  the 
purpose  of  power  management  control. 

Equivalent  circuit  battery  models  have  been  widely  used  for  the 
purpose  of  power  management  studies.  They  are  lumped  models 
with  relatively  small  number  of  parameters.  [42]  studies  a  number 
of  these  equivalent  circuit  models  for  Li-Ion  batteries.  As  the  battery 
SOC  in  the  current  research  is  maintained  within  a  narrow  range  of 


high  SOC  values,  the  simple  model  presented  in  Ref.  [42]  is  chosen. 
This  simple  model  is  a  SOC-dependent  voltage  source  in  series  with 
an  equivalent  resistance. 

The  above-mentioned  simple  model  is  used  only  for  the 
analytical  derivation  of  the  optimization  strategy.  However,  for  the 
system  simulation,  the  battery  model  of  the  Simulink  library  is 
used.  The  details  of  the  Simulink  battery  model  can  be  found  in  Ref. 
[43].  The  internal  resistance  of  a  Li-Ion  battery  bank  rated  at  80  V 
and  1000  Ah  is  0.8  mQ  which  is  calculated  through  the  user 
interface  of  the  Simulink  model.  Such  a  small  value  of  resistance  is 
resulted  from  a  high  number  of  parallel  branches  of  battery  cells  to 
provide  the  required  capacity. 

The  battery  loss  model,  based  on  its  equivalent  resistance  “Rb  \  is 
readily  integrated  in  (15)  by 

V  =  H  +  (16) 

Vbus 

where  Re ss  is  the  equivalent  resistance  of  the  energy  storage  system 
that  is  expressed  by 

n2 

^ess  =  ^bid+~ - Dy*^b  (^) 

The  above  simplified  model  of  the  energy  storage  system  is 
helpful  in  development  of  the  optimization  algorithm  that  will  be 
discussed  in  the  next  section. 
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Fig.  11.  The  power  loss  components  of  the  bidirectional  converter  in  boost  mode. 


Fig.  13.  The  converter  losses  based  on  the  equivalent  resistance. 
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3.  Optimization  algorithm 

The  presence  of  more  sources  of  power  causes  more  complexity 
in  power  management,  calling  for  an  effective  strategy  to  coordi¬ 
nate  the  sources  properly  and  to  make  an  appropriate  choice 
among  different  possible  alternatives.  Various  control  strategies 
can  be  devised  in  a  hybrid  power  system  so  as  to  achieve  desired 
performance  goals,  optimized  energy  efficiency,  reduced  emission, 
or  improved  reliability  and  safety.  This  section  deals  with  optimi¬ 
zation  of  energy  efficiency,  which  is  discussed  in  the  following. 

Fig.  14  compares  fuel  consumption  of  a  shipboard  dc  power 
system  with  that  of  a  conventional  shipboard  ac  power  system.  The 
solid  curves  depict  specific  fuel  consumption  for  the  dc  power 
system  of  Fig.  1,  wherein  “/<”  active  diesel  engines  are  controlled  at 
their  optimal  speeds.  This  set  of  active  generation  units  of  the  dc 
power  system  will  be  called  “dc  source”  hereafter  in  this  paper.  The 
dotted  curves  depict  specific  fuel  consumption  for  a  conventional 
ac  power  system,  wherein  the  same  set  of  diesel-generators  of  Fig.  1 
is  used,  and  “/<”  active  diesel  engines  operate  at  fixed  speed.  Besides 
the  comparison,  Fig.  14  also  implies  how  the  replacement  of  a  high- 
power  generation  unit  with  a  set  of  low-power  generation  units 
results  in  substantial  fuel  saving,  especially  under  light-load  con¬ 
ditions.  In  addition,  it  yields  a  quantitative  measure  to  determine 
the  optimal  number  of  active  generator  sets  based  on  the  system 
load.  However,  more  considerations  exist  with  regard  to  optimi¬ 
zation,  which  are  discussed  here  in  this  section. 

In  practice,  potential  large  transients  of  the  ship  electrical  load 
may  require  the  system  to  continuously  keep  a  redundant  active 
generator  set.  This  is  already  a  technical  requirement  for  several 
ship  types  including  OSVs,  which  leads  to  considerable  increase  in 
their  fuel  consumption.  In  case  an  energy  storage  device  is  inte¬ 
grated  in  the  system,  it  can  supply  the  transient  power  until  an 
engine  start-up  is  really  necessary.  This  would  provide  the  required 
redundancy  with  no  additional  fuel  consumption  24]. 

The  energy  storage  also  enables  the  engines  to  operate  at  their 
optimal  conditions.  It  is  apparent  in  Fig.  14  that  the  optimal  fuel 
consumption  is  achieved  when  the  engines  are  operated  at  their 
minimum  SFC  point.  However,  due  to  power  losses  of  the  energy 
storage,  the  optimal  operation  is  not  necessarily  achieved  by 
keeping  the  engine  operating  point  at  the  minimum  SFC.  A 
compromise  is  to  be  considered  between  the  minimum  engine  SFC 
and  the  maximum  efficiency  of  the  energy  storage  system.  Besides, 
existing  ripple  of  the  generator  load  causes  ripples  on  the  engine 
SFC  that  results  in  additional  fuel  consumption.  Therefore,  the 
power  ripple  of  the  generation  units  should  also  be  optimized.  This 
will  be  discussed  in  detail  later  in  this  section. 


Fig.  14.  SFC  curves  for  “k”  active  diesel  engines  in  a  dc  power  system  compared  to 
those  in  the  conventional  ac  power  system. 


According  to  the  above  discussion,  the  optimization  strategy  has 
two  main  tasks:  first,  to  calculate  the  optimal  value  for  the  average 
power  of  the  dc  source,  and  second,  to  find  the  optimal  value  for  the 
power  ripple  of  the  dc  source.  This  section  extends  the  work  of  [34] 
to  suit  the  multi-generator  shipboard  power  system  of  Fig.  1. 

The  power  of  the  dc  source  is  expressed  in  terms  of  an  average 
value  and  a  ripple  proportional  to  the  load  ripple. 

Ps(t)  =  Psa  +  a(PL(t)-PLa)  (0<a<l).  (18) 

where  Psa  and  PLa  are  the  averaged  values  of  the  dc  source  power 
and  the  load  respectively,  and  a  determines  the  fraction  of  the  load 
ripple  supplied  by  the  dc  source.  The  rest  of  the  load  ripple  is  to  be 
supplied  by  the  battery,  which  is  also  loaded  by  an  averaged  power 
of  Pu  -  Psa.  The  purpose  of  the  following  analysis  is  to  find  optimal 
values  for  Psa  and  a. 

3  A.  Charge— discharge  mode 

In  order  to  maintain  the  state  of  charge  (SOC)  of  the  battery 
within  a  certain  limit,  it  is  required  that  a  diesel  generator  starts 
when  SOC  reaches  the  lower  boundary,  and  shuts  down  as  soon  as 
SOC  exceeds  the  upper  limit.  Therefore,  two  operating  modes  can 
be  considered  for  the  dc  source,  which  are  addressed  here  by 
charge  mode  and  discharge  mode.  We  rewrite  (18)  for  the  two 
modes  of  discharge  and  charge  respectively. 

^s,k(0  =  Psa,k  +  «k(^(t)  “  Pla)  (°  <  ak  <  V-  H9) 

=  ^sa,/<+l  +  ^fc+1  (^l(0  —  ^La)  (9  <  ^k+1  <  1)  (20) 

where  k  is  the  number  of  active  engines  in  discharge  mode.  Fig.  15 
shows  an  overview  of  the  optimization  problem.  In  this  figure,  the 
system  is  supplying  a  load  with  an  average  of  Pl3  and  ripple  ampli¬ 
tude  of  PL.  Charging  and  discharging  the  energy  storage  periodically, 
we  can  move  the  operating  point  toward  lower  SFCs.  During  the 
charge  mode,  “/<  +  1”  active  engines  supply  the  load  while  charging 
the  battery  with  PCh,a-  This  moves  the  operating  point  from  point  A  to 
point  B,  reducing  the  SFC  value.  During  the  discharge  mode,  “k” 
engines  are  active  and  the  battery  is  discharged  with  PdiS,a  in  order  to 
assist  the  generators  to  supply  the  load.  The  operating  point,  in  this 
mode,  moves  to  point  B'.  In  order  to  find  the  optimal  work  points  B 
and  B',  the  optimal  average  value  of  the  dc  source  power  can  be 
calculated  for  both  charge  and  discharge  modes  i.e.  Psa,k  and  Psa,k+ 1- 

In  addition,  the  power  ripple  of  the  dc  source  may  also  cause 
additional  fuel  consumption.  It  is  seen  in  Fig.  15  that  although  the 
average  of  Psj<{t)  is  equal  to  Psaj0  the  average  of  SFC(PSfk(£))  is  greater 
than  SFC(Psa,k),  which  is  resulted  from  the  nonlinearity  of  SFC 
curves.  This  means  that  the  power  ripple  may  lead  to  additional 
fuel  consumption,  as  the  fuel  consumption  is  proportional  to  the 
average  of  SFC  over  time.  Thus,  the  power  ripple  of  the  dc  source 
should  be  optimized  as  well  as  its  average  power. 
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Based  on  the  foregoing  discussion,  the  set  of  variables  that  must 
be  optimized  consists  of  Psa<k,  Psa,k+i,  «k,  and  ak+\.  The  solution  of  the 
optimization  will  be  a  compromise  between  the  minimum  con¬ 
sumption  of  the  dc  source,  and  the  maximum  efficiency  of  the  energy 
storage  system.  The  number  which  is  the  number  of  active  en¬ 
gines  in  discharge  mode,  is  an  integer  between  0  and  3  for  the  system 
under  study.  It  is  determined  in  such  a  way  that  Pl3  lies  between  the 
minimum  SFC  points  of  the  two  neighboring  SFC  curves. 

Assuming  that  each  charge/discharge  cycle  lasts  for  Tc,  including 
a  charge  time  Tc h  and  a  discharge  time  Tdis,  we  obtain  the  total  fuel 
consumption  by  taking  integral  of  (1)  over  time  “t”. 


C(t)  —  Nsu|c0)/c+1Tch  +  J  Psjc+ i(t)dt 

+  bM  f  Ps  (t)dt  +  Cq  jcT^ is  +  ak  f  Ps  k(t)dt 
Tch  Tdis 

+  bk  j  Plk(t)dt  +  Csu  1.  (21) 

Tdis  J 


Considering  the  energy  balance  of  the  energy  storage  over  a 
charge/discharge  cycle,  we  can  obtain  the  different  time  periods  of 
Tc,  TCh,  and  TdiS  in  terms  of  the  other  variables  (Appendix  [see  5.b)], 
which  gives 


C(t)  =  t  j  C0  fc+1Ds  +  CokD'  +  aPLa  +  aDs^p-  [(Psafc+1  -  PLa)2 
l  ^bus 

+  Req(l  -  «k+ 1)2]  +  [(Psa,k  -  PU)2 

^bus 

+  PeqO  _  a/f)2]  +  t>k  i  lA>  (Psa.k+1  +  ait+lPeq) 

+  bkD '  (p|a +  ttkRgq)  + 

(24) 

where  Ds  =  1  —  D(  =  rch/Tcis  calculated  by 


Pu  -  Psa,fc  +  fr  [ (Psa,fc  -  Pta)2  +  %q(  1  -  «fc)2j 

_  bus  L _ J _ 

^sa,/<+l  —  Psajc  ~  vp  [(^sa,fc+l  —  ^sa,/<)  (^sa,/<+l  +  f’sa,^  —  ^Pla)  +  ^eq(^/c+l  —  ak)(afc+l  +  0Lk  —  2)j 


(25) 


where  Nsu  is  the  number  of  engine  start-ups  over  the  time  “t”.  An 
equivalent  value  for  fuel  consumption  at  each  engine  start-up  is 
also  considered  which  is  represented  by  Qu- 

The  equivalent  power  ripple  is  defined  by  means  of  an  inte¬ 
grative  function  Req  that  is  expressed  by 


The  fuel  consumption  represented  by  (24)  has  a  minimum  that 
satisfies 


*=0, 

dak  daM 


=  0, 


9C 


9  Pc 


=  o, 


9C 


sa  ,/< 


9P< 


=  0. 


sa,/c+l 


(26) 


Peq  — 


\f 


J  (PL(0  -  Pla)2dt. 


t 


(22) 


Taking  the  steps  provided  in  the  Appendix  [see  5.a)],  we  can 
derive  (23)  from  the  integral  terms  of  (21),  which  yields  the  fuel 
consumption  as 


C(t)  =  Nsu 


Q),/<+l^ch  +  ^0,/c^dis  +  aPLa(Pch  +  Tdis) 


Fig.  16  depicts  the  hourly  fuel  consumption  of  the  dc  power 
system  under  a  certain  load  in  terms  of  the  dc  source  loading 
conditions.  As  the  hourly  fuel  consumption  has  a  five  dimensional 
surface,  it  is  drawn  in  terms  of  only  one  pair  of  (a,  Psa)  at  each  time, 
assuming  the  other  pair  to  be  fixed  at(a*,P*a).  Therefore,  two 
separate  graphs  are  used  to  depict  the  fuel  consumption  against  the 
two  pairs.  The  optimal  values  of  Psa,k ,  Psa,k+ i,  a.k,  and  ak+\  can 
numerically  be  obtained  by  searching  for  the  minimum  fuel  con¬ 
sumption  in  the  analytical  relation  (24). 


+  aTch^  [(Psa.k+1  -  PLa)2  +  ^q(l  -  «fc+l)2] 

Vbus 

+  [(Psa,/<  -  PLa)2  +  ^q(l  -  «,<)2] 

Vbus 

+  bk+ 1  Tch  (^sa,/<+l  +  al+  l^eq)  +  bkTd[s  (P^ak  +  0^kReq) 

+  csu|. 

(23) 


3.2.  Continuous  mode 

The  obtained  optimal  loading  conditions  must  also  satisfy  the 
limits  on  the  engine  duty  cycle  “Ds”.  The  resulted  duty  cycle  must  be 
a  value  between  0  and  1 ;  otherwise  the  dc  source  operates  at  a 
continuous  mode.  In  case  the  resulted  duty  cycle  is  less  than  or 
close  to  zero,  the  dc  source  operates  with  “/<”  active  diesel  engines 
continuously,  and  if  the  duty  cycle  is  greater  than  or  close  to  one, 
the  dc  source  continuously  operates  with  “/<  +  1”  active  engines. 

In  such  conditions,  the  optimal  values  can  be  analytically 
attained  by 
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(a)  (b) 

Fig.  16.  Hourly  fuel  consumption  of  the  dc  source  under  a  certain  load,  a)  in  terms  of  loading  conditions  in  discharge  mode,  b)  in  terms  of  loading  conditions  in  charge  mode. 
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computational  intensity  in  search  for  optimal  ctk  and  au+v  The 

(27)  number  of  steps  in  search  for  PSSLtk  and  Psa,k+i  is  considered  20  in 
both  online  and  offline  processing.  Another  difference  between  the 
offline  algorithm  and  the  online  control  is  related  to  estimation  of 
the  load  average  “PLa”,  and  the  equivalent  power  ripple  “Req”.  As  the 

(28)  future  data  for  these  variables  are  not  available,  their  online  esti¬ 
mate  can  be  obtained  through  filtering  their  instantaneous  values 
by 


where  j 


/  k,  Ds  <  0 
\k  + 1,  DS>1 


In  the  following  section,  the  presented  algorithm  will  be  used  in 
an  online  optimization  strategy  for  an  offshore  support  vessel  at  its 
different  modes  of  operation. 


fUs)  - 

ReqM  =  r^(PL(s)-PLa(s))2. 


(29) 

(30) 


4.  Online  optimization  strategy 

The  optimization  algorithm  presented  in  the  previous  section 
should  be  tailored  to  suit  online  computations  and  control.  In  the 
theoretical  approach,  analytical  definitions  for  some  variables  are, 
as  such,  dependent  on  the  future  data  of  the  system,  for  example 
the  equivalent  power  ripple  “Req”  and  the  average  power  PLa.  As 
future  data  is  not  available  in  online  computations,  these  variables 
need  to  be  estimated.  In  addition,  the  online  processing  requires 
reduced  computational  intensity  to  perform  the  tasks  in  real  time. 
This  can  be  provided  by  reducing  the  step  sizes  in  online  processing 
when  searching  for  optimal  values  of  Psa>/<,  Psa,k+i,  and  ctk+ 1- 
These  estimations  and  simplifications  may  deteriorate  the  effec¬ 
tiveness  of  the  method.  In  order  to  evaluate  the  performance  of  the 
method  under  online  operation,  this  section  presents  the  simula¬ 
tion  results  of  the  online  control  in  comparison  with  those  of  the 
offline  theoretical  algorithm. 

An  online  control,  based  on  the  theoretical  algorithm  is  used  in 
simulation  of  an  offshore  support  vessel  over  a  number  of  opera¬ 
tional  modes.  The  parameters  of  the  system  are  reported  in  the 
Appendix.  The  power  demands  at  each  of  the  operating  modes  are 
specified  in  Table  1,  where  the  number  of  active  engines,  average 
power,  load  fluctuation,  and  duration  for  each  mode  are  given. 
These  data  are  deduced  from  the  experiences  of  ABB  and  DNV  re¬ 
ported  in  Refs.  [7,24]. 

For  each  mode  of  operation,  the  algorithm  outputs  the  optimal 
loading  conditions  of  the  dc  source,  i.e.  Psa>/<,  Psa,k+ i,  and  a^+i,  and 
also  the  control  mode  (charge-discharge  or  continuous).  The  al¬ 
gorithm  suggests  the  continuous  mode  when  the  operating  point  of 
the  dc  source  is  close  to  the  minimum  SFC;  otherwise  it  suggests 
the  charge— discharge  mode  as  the  optimal  mode.  While  the  offline 
algorithm  uses  a  step-size  of  0.01  in  search  for  optimal  au  and  au+b 
requiring  100  steps,  the  online  algorithm  uses  a  step-size  of  0.05, 
which  needs  20  steps.  This  provides  significant  reduction  in 


where  the  filter  time  constant  “T  should  properly  be  chosen;  A  too 
short  time  constant  causes  unnecessary  engine  start-stops,  when  a 
change  occurs  in  the  load  profile.  A  too  long  time  constant  causes 
that  the  dc  source  operates  under  non-optimal  conditions  for  an 
unnecessary  long  period  of  time  after  a  step  change  occurs  in  the 
system  load.  The  choice  of  the  time  constant  is  critical  considering 
that  the  load  variations  might  have  extremely  low  frequencies  in 
ship  applications,  e.g.  in  low  dynamic  positioning  mode.  In  order  to 
adapt  to  various  situations,  two  different  time  constants  have  been 
considered  in  the  simulations.  When  a  change  happens  in  the  load 
profile,  e.g.  changes  in  ship  operating  mode,  the  system  selects  the 
short  time  constant  throughout  the  transition  time  to  follow  the 
load  average  quickly.  In  the  normal  operation  at  each  mode,  the 
control  system  uses  the  larger  time  constant  that  results  in  more 
efficient  control.  The  mode  transition  is  detected  by  observing  a 
monotonic  change  in  the  load  average  for  some  seconds.  Repeated 
simulations  have  shown  that  suitable  values  for  the  short  and  long 
time  constants  are  300  s  and  1000  s  respectively.  Fig.  17  shows  the 
simulation  results  of  the  optimized  shipboard  dc  power  system  at 
three  modes  of  harbor,  high  dynamic  positioning,  and  transit  sup¬ 
ply.  It  is  seen  that  in  harbor  and  transit  supply  modes,  one  of  the 


Table  1 

Parameters  of  the  ship  operating  modes. 


Operating  mode 

#  of  active 
engines 

PLa  (%) 

Plp  (%) 

Duration  (h) 

Low  dynamic  positioning  (LDP) 

2 

33 

20 

7.5 

High  dynamic  positioning  (HDP) 

3 

53 

20 

4 

Anchor  handling  (AH) 

3 

50 

2 

1.5 

Harbor  (H) 

1 

17 

1 

1.5 

Bollard  pull  (BP) 

4 

66 

2 

0.5 

Transit  towing  (TT) 

4 

72 

1 

3 

Transit  supply  (TS) 

2 

33.5 

1 

6 

350 


B.  Zahedi  et  al.  /  Journal  of  Power  Sources  255  (2014)  341-354 


Harbor  Mode 


High  Dynamic  Positioning 


0  12  3  4 


Cont.: 

Ch./Dis.: 


12  3  4 

Time  (hr) 


Transit  Supply  Mode 


Fig.  17.  Simulation  results  for  three  operating  modes  of  the  shipboard  hybrid  power  system  with  online  optimization. 


diesel  engines  periodically  turns  on  and  off.  This  is  the  control 
mode  that  we  earlier  called  “charge-discharge  mode”.  The  battery 
state  of  charge  (SOC)  is  bounded  between  60  and  80%.  The  energy 
storage  system  periodically  charges  and  discharges  with  optimum 
average  values,  and  also  supplies  an  optimum  fraction  of  the  power 
ripple.  In  high  dynamic  positioning  mode,  in  which  the  load 
average  is  close  to  the  rated  conditions,  the  system  operates  at 
“continuous  mode”,  and  the  energy  storage  system  only  supplies  an 
optimum  fraction  of  the  power  ripple.  At  the  bottom  of  Fig.  17,  there 
is  a  flag  indicating  whether  the  control  mode  is  “charge-discharge” 
or  “continuous”. 

Simulation  results  over  an  operating  cycle  consisting  seven 
different  modes  are  presented  in  Fig.  18.  The  operating  modes  over 
the  cycle  are  indicated  on  top  of  the  figure.  The  upper  plot  shows 
the  instantaneous  and  average  load  power.  The  dual  time-constant 
characteristic  of  the  filter  is  particularly  clear  in  low  dynamic 
positioning  (LDP).  At  the  beginning  of  the  LDP  mode,  short  after  the 


mode  transition,  the  average  power  follows  the  instantaneous 
value,  as  the  smaller  time  constant  is  chosen.  Thereafter,  when  the 
transition  is  completely  over,  the  larger  time  constant  is  selected  in 
order  to  filter  the  extremely  low  frequency  variations.  This  proce¬ 
dure  is  carried  out  at  all  mode  transitions. 

During  harbor,  transit  towing,  and  transit  supply,  the  control 
system  operates  at  charge-discharge  mode,  as  periodical  changes 
in  the  number  of  active  engines  and  the  battery  SOC  are  seen. 
During  the  rest  of  ship  operating  modes,  the  control  system  oper¬ 
ates  at  continuous  mode,  and  the  battery  mainly  contributes  in 
supplying  the  power  ripple. 

Table  2  presents  the  optimized  fuel  consumption  in  comparison 
with  non-optimized  dc  and  ac  systems.  It  also  compares  fuel  con¬ 
sumption  results  for  the  online  and  the  offline  optimization 
methods  at  different  ship  operating  modes.  Regarding  specific  fuel 
consumption  (SFC),  two  modes  of  ATI  and  FIDP  have  average  power 
values  of  400  kW  which  are  in  the  proximity  of  the  minimum  SFC 


Fig.  18.  Simulation  results  of  the  optimized  shipboard  hybrid  power  system  over  an  operation  cycle  with  7  modes. 
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Table  2 

Fuel  consumption  based  on  the  algorithm  and  the  simulation,  compared  to  non-optimized  systems. 


Operating 

mode 

DC  system  with  ESS  (Optimized) 

DC  system  w/o  ESS 

AC  system 

Control 

mode 

#  of  active 
engines 

Aa.k+1  (%) 

P sa.fc  (%) 

ai<+i 

OL  k 

Copt  (gr/h) 
(algorithm) 

Copt  (gr/h) 
(simulation) 

Cdc  (gr/h) 

Cac  (gr/h) 

LDP 

Cont. 

1 

66.0 

N/A 

0.39 

N/A 

3.19E+04 

3.32E+04 

3.79E+04 

4.32E+04 

HDP 

Cont. 

2 

N/A 

79.0 

N/A 

0.49 

7.59E+04 

7.69E+04 

7.95E+04 

8.51  E+04 

AH 

Cont. 

2 

N/A 

75.0 

N/A 

0.49 

7.17E+04 

7.22E+04 

7.58E+04 

8.17E+04 

H 

Ch./dis. 

1/0 

62.8 

N/A 

0.36 

N/A 

9.83E+03 

1.16E+04 

1.35E+04 

1.69E+04 

BP 

Cont. 

3 

N/A 

75.3 

N/A 

0.78 

1.26E+05 

1.27E+05 

1.28E+05 

1.33E+05 

TT 

Ch./dis. 

4/3 

77.2 

94.6 

0.77 

0.72 

1.38E+05 

1.38E+05 

1.38E+05 

1.43E+05 

TS 

Ch./dis. 

2/1 

71.2 

100 

0.61 

0.44 

6.47E+04 

6.50E+04 

7.64E+04 

8.68E+04 

(see  Fig.  14).  One  may  argue  that  the  loading  points  of  this  study 
have  by  chance  fallen  in  inefficient  regions  so  that  the  resulted 
improvement  from  the  control  system  would  be  unrealistically 
high.  Although  the  load  profile  is  based  on  actual  data  reported  in 
Refs.  [7,24],  to  include  more  efficient  modes,  we  also  study  two 
modes  that  are  shifted  to  the  minimum  SFC  point  rated  at  800  kW. 
Let  call  these  modes  high  Bollard  Pull  (HBP)  and  High  Transit 
Towering  (HTT).  In  case  the  two  modes  BP  and  TT  are  replaced  with 
HBP  and  HTT  respectively,  4  out  of  7  operational  modes  are  in  the 
proximity  of  minimum  SFC  according  to  Fig.  14.  The  simulation 
results  for  these  two  modes  are  shown  in  Fig.  19.  It  is  seen  that 
under  this  loading  condition,  four  diesel  engines  operate  in 
continuous  mode,  and  the  battery  only  contributes  in  supplying  the 
power  ripple.  The  optimized  fuel  consumption  is  1.526e+5  (gr/h), 
which  is  exactly  the  same  as  that  of  the  dc  system  without  energy 
storage.  Therefore,  the  energy  storage  cannot  provide  any  fuel 
saving  for  the  dc  system  in  these  modes.  However,  the  fuel 


High  Transit  Towing  (HTT)  /  High  Bollard  Pull  (HBP) 
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Fig.  19.  Simulation  results  for  the  modes  HBP  and  HTT  of  the  dc  hybrid  power  system 
with  online  optimization. 


consumption  of  the  ac  system  is  1.567e+5  (gr/h)  under  the  same 
loading  condition.  This  means  2.6%  fuel  saving  is  still  provided  by 
the  optimized  system  compared  to  the  conventional  ac  system. 
These  fuel  saving  results  will  also  be  used  later  in  the  concluding 
evaluations  of  this  section. 

The  fuel  consumption  of  the  simulation  is  slightly  higher 
compared  to  the  algorithm,  which  is  caused  by  estimations  of  (29) 
and  (30)  as  well  as  the  step-size  reductions  in  search  for  the 
optimal  loading  conditions.  The  fuel  saving  achieved  by  the 
optimization  method  is  shown  in  percentage  in  Fig.  20.  The 
optimization  method  provides  more  fuel  saving  in  those  modes 
whose  power  demand  is  away  from  the  rated  conditions.  These 
modes  include  harbor,  transit  supply,  and  low  dynamic  posi¬ 
tioning  modes.  Due  to  approximations  made  in  the  online  opti¬ 
mization,  the  resulted  fuel  saving  of  the  simulation  is  less  than 
that  of  the  algorithm  by  around  1%  in  average.  The  difference 
between  the  performances  of  the  simulation  and  the  algorithm 
peaks  at  harbor  mode.  This  is  mainly  resulted  from  the  steep  SFC 
characteristic  over  the  low-power  range  which  the  harbor  mode 
belongs  to.  This  makes  the  harbor  mode  more  sensitive  to  the 
errors  resulted  from  the  approximations  of  the  online  processing. 
However,  due  to  the  low-power  characteristic  of  the  harbor  mode, 
it  has  a  minor  effect  on  the  overall  performance  of  the  online 
control. 

The  presented  results  in  Fig.  20  show  an  average  fuel  saving  of 
15.3%  for  online  control  simulation  when  the  conventional  ship¬ 
board  ac  power  system  is  replaced  by  a  dc  hybrid  power  system. 
This  percentage  is  calculated  based  on  weighted  averaging  of 
consumed  fuel  of  the  seven  different  operating  modes.  The  system 
simulation  over  the  24-h  cycle  results  in  even  a  slightly  higher  fuel 
saving  of  15.8%. 

Fig.  20  also  shows  that  integration  of  an  energy  storage  system 
into  the  dc  power  system  can  provide  a  fuel  saving  of  7.3%.  In  case 
the  modes  HBP  and  HTT  are  considered  instead  of  BP  and  TT 
respectively,  the  average  fuel  saving  would  be  14.6%  compared  to 
the  conventional  ac  power  system,  and  7.0%  compared  to  the  dc 
system  without  energy  storage.  It  should  be  noted  that  a  noticeable 
part  of  fuel  saving  is  resulted  from  optimizing  the  number  of  active 
generation  units.  As  addressed  in  Section  3,  the  conventional  OSVs 
are  required  to  continuously  keep  a  redundant  active  engine  during 
most  modes  of  operation  [see  Table  1].  However,  the  onboard  en¬ 
ergy  storage  of  a  hybrid  power  system,  while  optimizing  the 
loading  conditions  of  the  engines,  plays  the  role  of  a  redundant 
active  engine  that  enables  optimization  of  the  number  of  active 
engines  as  well.  Therefore,  the  achieved  fuel  saving  is  resulted 
partly  from  optimization  of  the  number  of  active  engines  and  partly 
from  optimization  of  loading  conditions  of  the  diesel-generators. 
The  study  results  demonstrate  a  significant  potential  regarding  fuel 
saving  in  offshore  support  vessels  through  transition  from  ac  to  dc 
power  systems  and  integration  of  energy  storage  devices.  The  fuel 
saving  potential  is  not  limited  to  the  studied  type  of  ships,  but  can 
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Fig.  20.  Resulted  fuel  saving  from  the  optimization  strategy,  a)  algorithm  results,  b)  online  simulation  results. 


be  also  considered  for  a  wide  range  of  marine  vessels  that  have 
frequent  load  variations  or  long-time  operation  under  non-rated 
conditions. 

5.  Conclusion 

In  this  paper,  detailed  efficiency  analysis  of  a  shipboard  dc 
hybrid  power  system  was  carried  out  in  order  to  evaluate  its  po¬ 
tential  fuel  saving.  It  was  shown  that  maximum  fuel  saving  con¬ 
ditions  for  such  a  system  required  the  generator  sets  to  operate  at 
optimal  loading  conditions  in  terms  of  average  power  and  power 
ripple.  A  numerical  optimization  algorithm  was  proposed  to  find 
optimal  loading  conditions,  using  the  presented  analysis.  In  addi¬ 
tion,  the  optimization  algorithm  determined  whether  the  gener¬ 
ator  sets  operated  at  continuous  mode  or  periodical  charge/ 
discharge  mode.  In  order  to  evaluate  the  actual  fuel  saving  pro¬ 
vided  by  the  method  when  processing  online,  the  algorithm  was 
used  in  simulation  of  an  Offshore  Support  Vessel  (OSV)  over  seven 
operating  modes.  Although  the  online  optimization  had  slightly 
less  efficiency  than  the  algorithm  by  around  1%,  it  provided  an 
average  fuel  saving  of  around  15%  compared  to  the  conventional  ac 
systems.  The  study  also  showed  that  the  shipboard  dc  system  with 
energy  storage  could  provide  a  fuel  saving  of  around  7%  compared 
to  the  dc  system  without  energy  storage.  The  study  results 
demonstrate  significant  fuel-saving  potentials  in  maritime  sector, 
particularly  in  the  category  of  marine  vessels  that  experience 
frequent  load  variations  or  long-time  operation  under  non-rated 
conditions. 
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Appendix 

1)  SFC  characteristics: 

l.a.  The  support  points  for  the  specific  fuel  consumption  of  a 
300  kW  diesel  engine:  Table  A.l 

1. b.  The  quadratic  function  coefficients  of  fuel  consumption  in 

terms  of  electrical  power:  Table  A.2 

2)  The  parameters  of  the  generation  units: 

2. a.  Diesel  engine:  6  cylinders,  30  rps,  Ky  =  1.5e6,  300  kW. 

2.b.  The  synchronous  generator:  300  kVA,  320/460  V,  2  pole 

pairs,  J  =  3.35  kg.m2,  rs  =  16.6  mQ,  =  5.245  mQ, 

Lifd  =  0.68  mH,  L\s  =  16.8  mH,  Lmd  =  5.85  mH, 

Lmq  =  5.05  mH,  r'kdi  =  0.1526  Q,  L'ikdi  =  3.404  mH, 

L'lkqi  =  3.404  mH,  r'kqi  =  40.57  mQ,  Cpe  =  3.4e-6,  cm  =  0.15. 

2.c.  The  rectifier:  VfW  =  1.5  V,  RD,on  =  1.2  mQ. 

2.d.  The  PI  controller  of  the  governor:  K P  =  0.0077,  Kj  =  0.0165. 

2. e.  Voltage  regulator-exciter:  I<P  =  0.4,  Kj  =  0.2,  Tfb  =  0.01  s, 

Iff  =  1  ms. 

3)  The  parameters  of  the  energy  storage  system: 

3. a.  The  low-voltage  side  of  the  dc/dc  converter:  VLV  =  100  V, 

L  =  516  pH,  Ci  =  500  pF,  n  =  3,  rL  =  3  mQ,  rCi  =  5  mQ, 
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rT  =  3  mQ,  RSi0n  =  0.6  mQ,  C0  =  50  nF,  Vfw  =  1.2  V, 
i^D.on  =  1.2  mQ. 

3.b.  The  high-voltage  side  of  the  dc/dc  converter:  Vhv  =  600  V, 
C2  =  267  pF,  rc2  =  0.34  Q,  Rs.on  =  2.9  mQ,  Co  =  2.3  nF, 
VW  =  0.8  V,  PD,on  =  0.05  mQ. 

3.c.  Li-Ion  battery  parameters:  Vfuii  =  93.1  V,  Vnom  =  80  V,  1  kAh, 
Rb  =  0.8  mQ. 

3.d.  The  ESS  loss  model  parameters:  % d  =  205  mQ,  Ress  =  0.3  Q, 
P10  =  100  W. 

4)  Loss  analysis  for  the  boost  mode: 


PSw,o 


=  J  (2Rs,on||2Rs,on)i?,  (0  <  t  <  DT) 

\2Rs,onPL,  (DT  <t<T) 


fsw\r  x/2  1 


vi 


hv 


rSw,sw  =  4(^1  C0V|M  =  2/swCb-  2  • 


Pd  —  2Vfw/bv  +  2R[),on4iv,rms 
=  2^3*  +  2RDon  f  Phv 

''hv  V^hvV  1  -  D 


Pl  =  rjlrms  =  rL 


nPhv 


(1  ~D)VhJ  ■ 


Pci  =  rCl  xl.rms  =  rCl 


PC2  =  rC2^C2,rms  =  rC2 


DP, 


hv 


(l-0)Vhv 


DVh 


V 1  2/swi  7 


PT  =  rT/f  =  rT 


nP 


hv 


(A.1) 


Psw.on  =  RsonllD  +  2%on^ (1  -  D)  =  ( 2  -  D)Ps,on^ •  (A.2) 


(A.3) 


(A.4) 


(A.5) 


(A.6) 


(A.7) 


(A.8) 


/  [*Wi(f)-PL(t)]dr-  J  P*s(t)dt 

Tch  Td  h 

+  /  [Ps,k(t)-PL(0]dt-  J  p^ss(t)dt 


=  0. 


(A.1 1 ) 


J  fWi(t)dt+  J  PsM(t)dt=  J  PL(t)dt+  J  PL(t)dt 

Tch  Tdis  Tch  rdis 

+  J  Pfhessm+  J  PfgsMt. 

TCh  Tdis 

(A.1 2) 

where  the  losses  of  the  energy  storage  system  are  obtained  by  using 
(16)  which  yields 


PfessiO  =^(PsM^)-PL(t))2. 

bus 

pfLd)  =  &(psm-pdt))2- 


(A.13) 


The  constant  term  of  the  ESS  losses  in  (16)  constitutes  a  power 
loss  of  100  W.  This  constant  value  is  considered  as  an  additional 
load  to  the  system  loading  so  that  only  the  second-order  term  will 
be  included  in  the  fuel  consumption  analysis. 

In  order  to  obtain  (23)  from  (21)  the  following  relations  have 
been  used.  Similar  relations  hold  for  PSii<. 

J  Pik+ i(t)dt  =  +  /  (fWi(t)  -*W+i)2df.  (A.14) 

TCh  Tc  h 

J  (Ps,k+l(t)-Psa,/<+l)2dt=  J  <1  CPL(t)  -  Pu)2dt 


(A.15) 


=  “fc+l^h^q- 


J  (Ps,k+l(t)-Psa,fc+l)(PL(t)-PLa)dt 
Tch 

=  J  ak+Mt)-Pla)2dt  =  aMTchR2eq. 


(A.16) 


5)  Analytical  derivation  of  the  fuel  consumption: 

5.  a  The  energy  balance  of  the  dc  bus  over  a  charge-discharge 
cycle  yields 


J  Pbus(t)dt  -  J  Pbus(t)df+  J  Pbus(t)dt  =  0.  (A.9) 

Tc  Trh  Ldis 

where 


^bus(0  — 


_  /  Ps,k+ 1  (t)  -  PL(t )  -  Pfess(t)  0  <  t  <  DSTC 
Ps±(C)-PL(C)-P(gs(C)  DSTC  <t<Tc 


®k+ 1  J  Ps, k+l  (f)df+  Qfc  J  Ps,k(t) df  —  oPla (Tch  +  Tdis) 

TCh  Tdis 

+  a7'chfe1[(PSa,k+l-Pla)2 

vbus 

+  ^eq(l  _2afc+1  •  ttk  t)] 

+  ardis^[(Psaik-Pu)2 

vbus 

+  R2eq(\  -2ak  + ufj]. 


(A.17) 


Replacing  (A.10)  in  (A.9),  we  will  have 


(A.10) 


5.  b  The  energy  balance  of  the  ESS  over  a  charge-discharge  cycle 
is  used  to  find  (25)  by  using  the  following  relations. 
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AEch  +  AEdis  =  0.  (A.18) 

AEch=  )  [Ps.k+](t)-PL(t)}c\t-  j  Pf%s(t)dt 
TCh  Tc  h 

=  rj  *W+1  -  Pu  -  It  [('W+I  -  Pi,)2 

l  Vbus 

+  ^q(l-«fc+l)2]j-  (A.19) 

A£dis=  /  [Ps,k(t)-Pi(t)]dt-  j  Pffss(t)dt 

^dis  Tdis 

=  Td  JpsaM  -  PLa  -  ^ [(Psajt  -  Pu)2  +  Rl q(l  -  «fc)2]  }. 

I  Vbus  J 

(A.20) 


Table  A.l 

Specific  fuel  consumption  against  mechanical  power. 


SFC  (gr/kWh) 

85  kW 

150  kW 

220  kW 

300  kW 

Fixed  speed 

242 

224 

204.3 

206 

Variable  speed 

219 

206.4 

197 

206 

Table  A.2 

Quadratic  function  coefficients  of  fuel  consumption  against  electrical  power. 

Coefficients 

Co  (gr/h) 

a  (gr/kWh) 

b  (gr/kWh  kW-1) 

Fixed  speed 

12761.7 

92.38 

0.235 

Variable  speed 

8488.1 

115.65 

0.202 
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